home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / slarrf.z / slarrf
Encoding:
Text File  |  2002-10-03  |  4.5 KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAARRRRRRRRFFFF((((3333SSSS))))                                                          SSSSLLLLAAAARRRRRRRRFFFF((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLARRF - the initial representation L D L^T and its cluster of close
  10.      eigenvalues (in a relative measure), W( IFIRST ), W( IFIRST+1 ), ..
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SLARRF( N, D, L, LD, LLD, IFIRST, ILAST, W, DPLUS, LPLUS,
  14.                         WORK, IWORK, INFO )
  15.  
  16.          INTEGER        IFIRST, ILAST, INFO, N
  17.  
  18.          INTEGER        IWORK( * )
  19.  
  20.          REAL           D( * ), DPLUS( * ), L( * ), LD( * ), LLD( * ), LPLUS(
  21.                         * ), W( * ), WORK( * )
  22.  
  23. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  24.      These routines are part of the SCSL Scientific Library and can be loaded
  25.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  26.      directs the linker to use the multi-processor version of the library.
  27.  
  28.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  29.      4 bytes (32 bits). Another version of SCSL is available in which integers
  30.      are 8 bytes (64 bits).  This version allows the user access to larger
  31.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  32.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  33.      only one of the two versions; 4-byte integer and 8-byte integer library
  34.      calls cannot be mixed.
  35.  
  36. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  37.      Given the initial representation L D L^T and its cluster of close
  38.      eigenvalues (in a relative measure), W( IFIRST ), W( IFIRST+1 ), ... W(
  39.      ILAST ), SLARRF finds a new relatively robust representation L D L^T -
  40.      SIGMA I = L(+) D(+) L(+)^T such that at least one of the eigenvalues of
  41.      L(+) D(+) L(+)^T is relatively isolated.
  42.  
  43.  
  44. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  45.      N       (input) INTEGER
  46.              The order of the matrix.
  47.  
  48.      D       (input) REAL array, dimension (N)
  49.              The n diagonal elements of the diagonal matrix D.
  50.  
  51.      L       (input) REAL array, dimension (N-1)
  52.              The (n-1) subdiagonal elements of the unit bidiagonal matrix L.
  53.  
  54.      LD      (input) REAL array, dimension (N-1)
  55.              The n-1 elements L(i)*D(i).
  56.  
  57.      LLD     (input) REAL array, dimension (N-1)
  58.              The n-1 elements L(i)*L(i)*D(i).
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAARRRRRRRRFFFF((((3333SSSS))))                                                          SSSSLLLLAAAARRRRRRRRFFFF((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      IFIRST  (input) INTEGER
  75.              The index of the first eigenvalue in the cluster.
  76.  
  77.      ILAST   (input) INTEGER
  78.              The index of the last eigenvalue in the cluster.
  79.  
  80.      W       (input/output) REAL array, dimension (N)
  81.              On input, the eigenvalues of L D L^T in ascending order.  W(
  82.              IFIRST ) through W( ILAST ) form the cluster of relatively close
  83.              eigenalues.  On output, W( IFIRST ) thru' W( ILAST ) are
  84.              estimates of the corresponding eigenvalues of L(+) D(+) L(+)^T.
  85.  
  86.      SIGMA   (input) REAL
  87.              The shift used to form L(+) D(+) L(+)^T.
  88.  
  89.      DPLUS   (output) REAL array, dimension (N)
  90.              The n diagonal elements of the diagonal matrix D(+).
  91.  
  92.      LPLUS   (output) REAL array, dimension (N)
  93.              The first (n-1) elements of LPLUS contain the subdiagonal
  94.              elements of the unit bidiagonal matrix L(+). LPLUS( N ) is set to
  95.              SIGMA.
  96.  
  97.      WORK    (input) REAL array, dimension (???)
  98.              Workspace.
  99.  
  100.      IWORK   (input) INTEGER array, dimension (???)
  101.              Workspace.
  102.  
  103.      INFO    (output) INTEGER
  104.              = 0:  successful exit
  105.              > 0:  if INFO = 1, IFIRST and ILAST had illegal values
  106.  
  107. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  108.      Based on contributions by
  109.         Inderjit Dhillon, IBM Almaden, USA
  110.         Osni Marques, LBNL/NERSC, USA
  111.  
  112.  
  113. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  114.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  115.  
  116.      This man page is available only online.
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.